Non-equilibrium effects in transport through quantum dots 
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The role of non-equilibrium effects in the conductance through quantum dots is investigated. As- 
sociated with single-electron tunneling are shake-up processes and the formation of excitonic-like 
resonances. They change qualitatively the low temperature properties of the system. We analyze by 
quantum Monte Carlo methods the renormalization of the effective capacitance and the gate- voltage 
dependent conductance. Experimental relevance is discussed. 



PACS numbers: 73.40Gk, 73.23.Hk 
I. INTRODUCTION. 

Quantum dots are paradigms to study the transition 
from macroscopic to microscopic physics. At present, 
the role of single-electron charging is well understood 
I]. Processes which otherwise are found in solids at the 
single-atom level, such as the Kondo effect, are being cur- 
rently investigated . Other atomic features, like the ex- 
istence of high spin ground states have also been observed 
||. The existence of many internal degrees of freedom 
within the dot leads to a variety of effects reminiscent of 
those found in large molecules. Shake-up processes, asso- 
ciated with the rearrangement of many electronic levels 
upon the addition of one electron, have been reported 



In the present work, we will show that internal excita- 
tions of the dot lead to non-equilibrium effects which can 
substantially modify the transport properties. In suffi- 
ciently small dots, the addition of a single electron may 
cause significant charge rearrangements 0, and the re- 
sulting change in the electrostatic potential of the dot 
modifies the electronic level structure. In the limit when 
the level separation is much smaller than other rele- 
vant scales, this process leads to an "orthogonality catas- 
trophe" ||, first discussed in relation with the sudden 
switching of a local potential in a bulk metal. In addi- 
tion, an electron tunneling event changing the charge of 
the dot is associated with a charge depletion in the leads 
or in neighboring dots. The attraction between the elec- 
tron in the dot and the induced positive charge leads to 
the formation of an excitonic resonance, similar to the 
well known excitonic effects in X-ray absorption [p|JlO|] . 

The relevance of these effects for tunneling processes in 
mesoscopic systems was first discussed in refs. [ pl]Jr^ ] (see 
also jn|). The existence of excitonic-like resonances has 
indeed been reported in mesoscopic systems [ [l4|| . The 
formation of bonding resonances has been observed in 
double well systems Non-equilibrium effects like 

those studied here have also been discussed in Ref. fl^ j 
where transport through a tunnel junction via localized 
levels due to impurities was analyzed, and have been ob- 
served experimentally in Ref. |l7| |. A different way for 
modifying the conductance of quantum dots through the 



formation of excitonic states has been proposed in [JL8J. 
In that case the exciton is a real bound state, while the 
excitonic mechanism discussed here is a dynamical pro- 
cess. In some devices, the charge rearrangement may 
take place far from the tunneling region. In this case 
no excitons are formed but the orthogonality catastro- 
phe persists. This process has been discussed in relation 
to the measurement of the charge in a quantum dot by 
the current through a neighboring point contact |Q . 

We will study the simplest deviations from the stan- 
dard Coulomb blockade regime. Our analysis is valid for 
quantum dots where the spacing between electronic lev- 
els, Ae, is much smaller than the charging energy, Ec, 
or the temperature T . The non-equilibrium effects dis- 
cussed here will be observable if, in addition, the number 
of electrons in the dot is not too large so that changes 
in the charge state lead to non-uniform redistributions of 
the charge. Thus, we will consider an intermediate situ- 
ation between the Kondo regime and orthodox Coulomb 
blockade (see below for estimates). The processes that 
we consider will be present in double-dot devices, but, 
for definiteness, we study here a single dot. 

The main features discussed above can be described 
by a generalization of the dissipative quantum rotor 
model p^l , which has been studied widely in connection 
with conventional Coulomb blockade processes [^0|. We 
present a detailed numerical analysis of the generalized 
model, along similar lines as previous work by some of us 
on the dissipative quantum rotor 

The paper is organized as follows: In the next section, 
we show how to estimate the parameters which charac- 
terize non-equilibrium effects. The model is reviewed in 
section III, with emphasis on details needed for the sub- 
sequent calculations. The numerical method is presented 
in section IV. In Sees. V and VI, we give results for the 
renormalized capacitance and conductance. In Sec. VII 
we discuss some possible experimental evidence of the 
effects studied here. We close with some conclusions. 



II. NON-EQUILIBRIUM EFFECTS. 
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A. Inhomogeneous charge redistribution. 



B. Effective tunneling density of states. 



The standard Coulomb blockade model assumes that, 
upon a change in the charge state of the dot, the elec- 
tronic levels within a quantum dot are rigidly shifted by 
the charging energy, Eq — e 2 /2C with C the capacitance 
of the dot . Deviations from this assumption have been 
studied by means of an expansion in terms of g -1 , the 
inverse dimensionless conductance, g ~ k-pl, where fcp is 
the Fermi wave- vector, and I is the mean free path H| . It 
is also assumed that kp is small compared to the inverse 
Thomas-Fermi screening length, /ctf = \/^e 2 N(eF)/eo. 
To lowest order beyond standard Coulomb blockade ef- 
fects, the change in the charge state of the dot leads to 
an inhomogeneous potential, and induces a term in the 
Hamiltonian, which can be written as (2^J|j|: 

Hint = {Q~ Qoffset) / ^{v)U{v)4,{v)d d V . (1) 

Here Qoffset denotes offset charges in the environment, 
the operator Q measures the total electronic charge in 
the dot, and if>' (?) creates an electron at position r. The 
potential U(r) modifies the constant shift of the energy 
levels of the dot assumed in the standard Coulomb block- 
ade model. It appears due to the restricted geometry of 
the dot. After a charge tunnels into the dot, a pile-up of 
electrons at the surface of the dot is induced. As a result 
there is a net attraction of electrons towards the surface, 
besides a constant shift given by e 2 /2C. In general, the 
potential U in eq. ([!]) can be obtained from the Hartree 
approximation (in the Thomas-Fermi limit) for dots and 
leads of arbitrary shape. For a spherical dot of radius R 
the potential U(r) has the simple form |22]: 



U(r) = - 



£ A;tf-R 2 

and, for a two-dimensional circular dot: 
U(r) = - 



K 



2e k T FRv r R T 



K 



(2) 



(3) 



K is a constant which ensures (U(r)) = 0. Non- 
equilibrium effects arise because the potential in eq. (|l|) is 
time dependent, as it changes upon the addition of elec- 
trons to the dot. Hamiltonians with terms such as eq. (|l|) 
were first discussed in Q (see next section). Note that 
the potential is localized in the surface region, where the 
tunneling electron is supposed to land. Finally, the po- 
tential is attractive, leading to the localization of the new 
electron near the surface, giving rise to excitonic effects. 

Other non-equilibrium effects can arise if the charge of 
the dot induces inhomogeneous potentials in other metal- 
lic regions of the device. In this case, the only effect ex- 
pected is the orthogonality catastrophe (see below), due 
to the shake-up of the electrons both in the dot and in 
the other regions. We take this possibility into account 
in the analysis in the following sections. 



In the absence of non-equilibrium effects, the conduc- 
tance of a junction between the dot and the leads is 



2f_ 

h 



J2 NXlcad(£F)^,dot(£ F ) 



(4) 



cha 



where the summation is over the channels, is the hop- 
ping matrix element through channel i, N(Ep) is the 
density of states at the Fermi level, and we use the stan- 
dard theory of tunneling in the weak transmission limit 
p4j . Eq. (H|) implicitly assumes a constant density of 
states, as appropriate for a metallic contact. 

The non-equilibrium effects to be considered can be 
taken into account through a modification of the effective 
tunneling density of states jl2],[25| . In this case the elec- 
tron propagators in eq. ^ are the non-equilibrium ones, 
in an analogous way to the modifications required in the 
study of X-ray absorption spectra of core levels in metals 
H , or tunneling between Luttinger liquids p6|-p8|| . As in 
those problems, we can distinguish two cases: 

i) The analogue of the X-ray absorption process: The 
charging of the dot leads to an effective potential which 
modifies the electronic levels. At the same time, an elec- 
tronic state localized in a region within the range of the 
potential is filled. The interaction between the electron 
in this state and the induced potential must be taken 
into account (the excitonic effects, in the language of the 
Mahan-Nozieres-de Dominicis theory). 

ii) The analogue of X-ray photoemission: The charging 
process leads to a potential which modifies the electronic 
levels. The tunneling electron appears in a region outside 
the range of this potential. Only the orthogonality effect 
caused by the potential needs to be included. 

Taking into account the distinctions between these two 
possibilities, the effective (non-equilibrium) density of 
states in the lead and the dot becomes (omitting the 
channel index, i): 



D e s(uj) 



OC \<jj\ 



(5) 



with e given by 

3=1,2 Z 7T 



Si 



e = < 



^3=1,2 ( X ) 



(excitonic 
resonance) 

(orthogonality 
catastrophe) 



(6) 



Here Sj is the phase shift induced by the new electrostatic 
potential in the lead states (j = 1) or in the dot states 
(J = 2). The exponent is positive, e > 0, if excitonic 
effects prevail, while e < if the leading process is the 
orthogonality catastrophe. 
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We can get an accurate estimate for e in the simple 
cases of a spherical or circular quantum dot decoupled 
from other metallic regions discussed in . We assume 
that tunneling takes place through a single channel, and 



the contact is of linear dimensions 



In Born ap- 



proximation the effective phase shift becomes: 



N(e F ) [ U{r)d d r : 
Jn 



e lpkf^ spherical dot 
S N ft£r circular dot 



(7) 



where f2 is the region where tunneling processes to the 
leads are non negligible, typically of dimensions compa- 
rable to kp . Note that, for a very elongated dot (d=l), 
the phase shift will not depend on its linear size. As men- 
tioned earlier, the leads can modify significantly these es- 
timates. The tunneling electron can be attracted to the 
image potential that it induces, enhancing the excitonic 
effects (e > 0). On the other hand, shake-up processes in 
metallic regions decoupled from the tunneling processes 
will increase the orthogonality catastrophe, without con- 
tributing to the formation of the excitonic resonance at 
the Fermi energy. 



III. THE MODEL. 

The shake-up processes mentioned in the preceding 
section are described by the Hamiltonian [113] : 



H = Hq + Hr, + Hl + Ti-T + H\. 



T~(-Q 



(Q-Q 



offset ) 



2C 



Hi = e k,iCk,i c k,i , « = L,R 
fc 

k,k' 



fc.fc' 



(8) 



Here [<fi, Q] — ie. The Hamiltonian separates the junc- 
tion degrees of freedom into a collective mode, the charge 
Q, and the electron degrees of freedom of the electrodes 
and the dot. This separation is standard in analyzing 
electron liquids, where collective charge oscillations (the 
plasmons) are treated separately from the low-energy 
electron-hole excitations. In our case, this implies that 
only those states with energies lower than the charging 
energy are to be included in Hi, Ht and Hint in eq.(||). 
Higher electronic states contribute to the dynamics of 
the charge, described by Hq. The Hamiltonian, eq.(||), 
suffices to describe transport processes at voltages and 
temperatures smaller than the charging energy. 



In the following, we will express the offset charge 
Ooffset) introduced in eq. ([[]), by the dimensionless pa- 
rameter n e = Q ffset/e- By Vk,k' we denote the ma- 
trix elements of U(r) in the basis of the cigenfunctions 
near the Fermi level. We allowed that inhomogeneous 
potentials can be generated on both sides of the junc- 
tion. We assume that tunneling can take place through 
several channels (index channel has been omitted). The 
transmission through each channel should be small for 
perturbation theory to apply. 

The electrical relaxation associated with the tunneling 
process takes place in two stages. In the first, the tun- 
neling electron is screened by the excitation of plasmons, 
forming the screened Coulomb potential. The time scale 
for this process is of the order of the inverse plasma fre- 
quency. Next, the screened Coulomb potential excites 
electron-hole pairs. As the electrons at the Fermi level 
have a much longer response time, they feel this change 
as a sudden and local perturbation. 

We will restrict ourselves to the regime where the level 
spacing is small, Ae <C T,E C . Using standard techniques 
[^,^9|, we can integrate out the electron- hole pairs and 
describe the system in terms of the phase </> and charge Q 
only. This procedure leads to retarded interactions which 
are long-range in time, as the electron-hole pairs have a 
continuous spectrum down to zero energy. It is best to 
describe the resulting model within a path-integral for- 
malism. Because of the non-equilibrium effects the effec- 
tive action is a generalization of that derived for tunnel 
junctions |1|9| 



S\d>] = 



dr- 



4£ c 



a 



13 



13 



dr I dr'E^ I - 



( 1 - cos[0(r) - <j>(r')} 
sin 2 - £ [7r(r-T')//3] 

(9) 

It describes the low energy processes below an upper cut- 
off of order of the unscreened charging energy, Ec- The 
parameter a cx t 2 Nu{e-p)NL{^F) is a measure of the high 
temperature conductance, go, in units of e 2 /h: 

9o 



47r 2 (e 2 //i) 



(10) 



Note that the definition of the action, eq. (^), does not 
allow us to study temperatures much higher than Eq- 
The kernel which describes the retarded interaction is 
given by the effective tunneling density of states, eq. (||). 
The value of e is the anomalous exponent in the tunneling 
density of states, given in eq. . 

The action ([)]) has been studied extensively for e = 
[20 2l], |30| , |3l|| , describing charging effects in the single- 
electron transistor in the usual limit where the electrodes 
are assumed to be in equilibrium. If e > 0, the model 
has a non-trivial phase transition |32|-|3^]. In this case, 
for a > OL CT it ~ 2/(7r 2 e), the system develops long-range 
order when T = l/j3 — ► 0, leading to phase coherence 
and a diverging conductance. 
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IV. COMPUTATIONAL METHOD 

For a given offset charge n e , the grand partition func- 
tion can be written in terms of the phase (f>, as a path 
integral [M: 



Z(n e 



exp(27nmn e ) / T>§ exp (— S[(j)]) , (11) 



where m is the winding number of 4>, and the paths 4>(t) 
satisfy in sector m the boundary condition </>(/?) = (f)(0) + 

The effective action and partition function can be 
rewritten in terms of the phase fluctuations 9(t) = 
4>(t) — 2-KrnT I with boundary condition 0((3) = 0(0), in 
the form 



Z— exp(27r«mrt e ) I m (a, e, (3) 



(12) 



The coefficients I m (a,e,(3) = f T>9 exp (—S m [8]) are 
to be evaluated with the effective action S m [6(T)] = 
S[9(t) + 27rmr//3]. They depend on the winding number 
m, the temperature, and the dimcnsionless parameters a 
and e, but are independent of the offset charge n e . This 
means that the problem reduces, from a computational 
point of view, to the calculation of the relative values 
of I m (a,e, /?), which can be obtained from Monte Carlo 
(MC) simulations apart from an overall normalization 
constant | 2l]j3l| ]. The partition function is even and pe- 
riodic with respect to n e , Z(n e ) = Z(—n e ) = Z(n e + 1), 
and therefore one can restrict the analysis to the range 

< n e < 0.5. 

The MC simulations have been carried out by the 
usual discretization of the quantum paths into N (Trot- 
ter number) imaginary-time slices fl3q| . In order to keep 
roughly the same precision in the calculated quantities, 
as the temperature is reduced, the number of time-slices 
N has to increase as 1/T. We have found that a value 
N = 4(3 Ec is sufficient to reach convergence of I m . 
Therefore, the imaginary-time slice employed in the dis- 
cretization of the paths is At = (3/N = l/(42?c). 

When discretizing the paths 4>(t) into N points, it is 
important to treat correctly the |r — t'| — > divergence 
that appears in the tunneling term St [<j>] of the effective 
action [second term on the r.h.s. of Eq. M)]. This di- 
vergence can be handled as follows. In the discretization 
procedure, the double integral in St[4>] translates into a 
sum extended to iV 2 two-dimensional plaquettes, each 
one with area (At) 2 . The above-mentioned divergence 
appears in the TV "diagonal" terms (t = r') and can 
be dealt with by approximating the integrand close to 
t = t' by E c \t - r'| e (d0/dr) 2 /2. Thus, by integrating 
this expression over the "diagonal" plaquette (i, i), with 

1 < i < N, one finds that its contribution to St[<f>] is 
given by 



AS t (r i ,T i ) = 2Eh— [ — ) \f 



At 

2 



c+2 



(13) 



which is regular for f / — 1. The error introduced by 
this replacement in the discretization procedure is of the 
same order as that introduced by the usual discretization 
of the "non-diagonal" terms. We have checked that the 
results of our Monte Carlo simulations obtained by using 
this procedure converge with the Trotter number N. 



has been sampled 
36| for temperatures 



The partition function in Eq.(]12|) 
by the classical Metropolis method 
down to k^T — Ec /200. A simulation run proceeds via 
successive MC steps. In each step, all path-coordinates 
(imaginary-time slices) are updated. For each set of pa- 
rameters (a, e, T), the maximum distance allowed for ran- 
dom moves was fixed in order to obtain an acceptance ra- 
tio of about 50%. Then, we chose a starting configuration 
for the MC runs after system equilibration during about 
3 x 10 MC steps. Finally, ensemble-averaged values for 
the quantities of interest were calculated from samples of 
~ 1 x 10 5 quantum paths. More details on this kind of 
MC simulations can be found elsewhere [Ell . 



V. RENORMALIZATION OF THE 
CAPACITANCE. 



We will study the capacitance renormalization for tun- 
neling conductance a > by calculating the effective 
charging energy Eq(T) = e 2 /2C*(T), which can be ob- 
tained as a second derivative of the free energy F = 
— fceTlnZ: 



= 2 m 



(14) 



n e =0 



At high temperatures, the free energy F(n e ) depends 
weakly on n e , and the curvature [i.e., E^(T)] approaches 
zero. At low T, and for weak tunneling (a <^ 1), it 
coincides with the usual charging energy Ec- By using 
Eqs. ([ll]) and (|l4| ) this renormalized charging energy can 
readily be expressed as 



E* C (T) = 2^ 2 fc B T(m 2 ) 



(15) 



where 



, = o is the second moment of the coefficients 



The correlation function in imaginary time G(t) , that 
will be used below to calculate the conductance, can 
be calculated from the MC simulations as G(r) = 
(cos[i/)(t) — 0(0)]). This means in our context: 
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FIG. 1. Effective charging energy of the single electron 
transistor for different values of e and a = 0.25, as a func- 
tion of temperature. 




FIG. 2. Effective charging energy of the single electron 
transistor for different values of e, as a function of a for 
T/E c = 0.02. 



1 

G(r) = — exp(27rzmn e ) x 



m— — oo 



£><?() e - SM cos[<£(t) - </>(0)] 



££(a,e) 



a C rit(e) 



(17) 



In this expression we use a cr i t (e) as the only adjustable 
(16) parameter. The results are shown in fig. (^). Note that 
the same equation ( |i~7| ) can be used to fit the results for 
e < 0, if one uses a negative value for a cr it- There is no 
phase transition, however, for e < 0. 



A number of features, directly related to the free en- 
ergy of the model given by action (^J) , are reasonably well 
understood for e = 0. The effective charge induced by an 
arbitrary offset charge, n e , has been extensively analyzed. 
A number of analytical schemes give consistent results in 
the weak coupling (a « 1) regime p7| , p8| . These calcu- 
lations have been extended to the strong coupling limit 
by numerical methods |21|. 



In fig. ([!]), we present results for the effective charging 
energy Eq as a function of temperature, and for different 
values of e. The value of E* c is enhanced for e < 0, where 
the orthogonality catastrophe dominates the physics. A 
positive e reduces the effective charging energy, and, be- 
yond some critical value (see discussion below), scales 
towards zero as the temperature is decreased, showing 
non-monotonic behavior. The same trend can be appre- 
ciated in fig. (||), where the effective charging energy at 
low temperatures is plotted as a function of a. Renor- 
malization group arguments |32|,[33| show that E^(a, e) 
should go to zero for a > a cr u(e). We have checked the 
consistency of this prediction with the numerical results 
by fitting the values of E^ia) at low temperatures by 
the expression expected from the scaling analysis near 
the transition 1331: 




FIG. 3. Critical line determined by the RG calculation: 
Oicrit = 2/(7r 2 e) and values calculated fitting the numerical 
data for E^{a, e) to the expression (|l7|). 

The results show that the present calculations are very 
accurate even for relatively large values of a, where Eq 
converges at very low temperatures. 
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VI. EVALUATION OF THE CONDUCTANCE. 



while, on the other hand: 



A. e = 

The conductance of the single-electron transistor is no- 
toriously more difficult to calculate than standard ther- 
modynamic averages. It cannot be derived in a simple 
fashion from the partition function, and requires the ana- 
lytical continuation of the response functions from imagi- 
nary to real times or frequencies | j39| . Hence, there are no 
comprehensive results valid for the whole range of values 
of a, T/Ec and n e . For e = 0, the conductance g(T) can 
be written as |4ll-|4T 



g(T) = 2g p 



duj luS(lo) 
2tt eP" - 1 



(18) 



where go is the normal state conductance, and S(u>) is 
related to the correlation function in imaginary time 



ir \/3 



(25) 



If uiq = 0, both g and goG((3/2) have the same limit, as 
T — > 0. When uj ^ the leading term goes as T 2 (cotun- 
neling) in both cases, with prefactors equal to 2A04A/ir 
and AA/tt, respectively. 

From the above discussion of the relation between the 
high- and low-temperature behavior of g(T) and G((3/2), 
we find that the interpolation formula 



g(T) « g G 



(26) 



G(t) 



/g^(-r) e ~i<t>(0)\ 



1 i-oo e (p - T )w + gT« 
— / rfw — : A(w) 



2tt 



;/3w _ ^ 



and 



= (1 - e-P")S(u) 



(19) 



(20) 



In the previous expressions, the charging energy is the 
natural cutoff for the energy integrals. 
_At high temperatures, 0Ec ~ 1, we can expand in eq. 



1 fa 2 

(3 24 



7^/ 
5760 



so that: 
g(T) fa ff0 



7/3 4 
5760 



,(0u>)/2 



A(co) 



e l3u _ 1 

(21) 



G w 



(22) 



At low temperatures, (3Ec 3> 1, the conductance is 
dominated by the low energy behavior of A{oj) or, alter- 
natively, S(u>). To lowest order, we expect an expansion 
of the form: 



S(lj) pa 2tt5{u - uj q ) + A\w\ + 



(23) 



where u>o is an energy of the order of the renormalized 
charging energy (or zero at resonance, n e pa 1/2), and A 
is a constant which describes cotunneling processes ||4[ . 
Inserting this expression in eq. (|l8"|), we obtain: 



g-(T) pa ff0 - 



2(3uj 



.90 



A 



da; 

e»-l 



(24) 



should give a reasonable approximation over the entire 
range of parameters (note that the above discussion is 
independent of the values of a and n e ). Eq. ( p6| ) is 
consistent with the the main physical features expected 
both in the high and low temperature limits, at and away 
from resonance. The advantage of using G(/3/2) is that it 
can be computed, to a high degree of accuracy, by stan- 
dard Monte Carlo techniques, as it does not require to 
continue the results to real times. A similar approxima- 
tion, used to avoid inaccurate analytical continuations 
has been applied for bulk systems in Ref. pa]. 



We show the adequacy of the approximation, eq. (|2£ 
by plotting the conductances estimated in this way, as a 
function of the bias charge n e , in fig. (0). 
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FIG. 4. Conductance of the single-electron transistor 
(e = 0), as a function of the dimensionless bias charge n e , 
for several values of the coupling a, and two different temper- 
atures. 
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FIG. 5. Maximum and minimum values of the conductance 
of the single-electron transistor (e = 0) , as a function of tem- 
perature, for different values of the coupling a. 

The minimum (n e — 0) and maximum conductances 
(n e = 1/2) for different values of a and temperatures are 
shown in fig. (||). 



B. e/0 

We now extend the previous approximation to the con- 
ductance, eq. (U6\j to the case f / 0. The main modifi- 
cation in eq. (O) is that a factor lo/{\ — er^ u ~) within 
the integral has to be replaced by the effective tunnel- 
ing density of states, -D c ff, given, at zero temperature, 
by eq. (||). At finite temperatures, the corresponding ex- 
pression is approximately 



Deff(w) cx max[T(T/E c )- e ,to(Lo/E c y 



(27) 



The relevant range in the integrand in eq. ( |l8| ) is from 
oj = to u> « T. In the following, we will factor the e 
dependent part of the effective density of states, and we 
write the generalization of eq. (^|) to finite temperatures 
as: 



1 - e-@ u 



(28) 



Finally, when inserting this expression into eq. (^), we 
make the approximation: 



D res (e,u) 



T 



(29) 



With this approximation, we can perform the same anal- 
ysis in the high and low temperature regimes as before, 
to obtain the interpolation formula: 



9(T) « go 



T 



G 



(30) 



This expression includes again the relevant physical pro- 
cesses at high and low temperatures. 



Results for the maximum and minimum values of the 
conductances, for different values of e, are presented in 
fig. In the non phase-coherent regime, at very low 
temperatures, T <C Eq, the conductance away from res- 
onance should vary as g cx T 2 ~ 2e . Exactly at resonance, 
n e = 1/2, the conductance diverges as T~ e . The most 
interesting result is the divergence of the conductance, 
at low temperatures, for e = 0.5, where the excitonic ef- 
fects are strong enough to drive the system to the phase- 
coherent phase. 
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FIG. 6. Maximum and minimum values of the conductance 
of the single-electron transistor, for different values of e. 



The full conductance, as a function of n e , is shown in 
fig. (0), for a — 0.25 and e = 0.5. As mentioned above, 
for these parameters the system is already in the phase 
coherent regime. The conductance behaves in a way sim- 
ilar to that in the usual case (e = 0), and a peaked 
structure develops. The absolute magnitude, however, 
increases as the temperature is lowered. Note that the 
effective charging energy is finite as T — > (see fig. (0)). 
It is interesting to note that, in this phase with complete 
suppression of Coulomb blockade effects at low tempera- 
tures (high values of e and high conductances), the peak 
structure appears only for an intermediate range of tem- 
peratures, and it is washed out at very low temperatures. 
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FIG. 7. Conductance, as a function of n e and temperature, 
of a single-electron transistor with a — 0.25 and e = 0.5. 

In fig. (||) we show the conductance as a function of e 
for a fixed temperature. It is evident the increase of the 
conductance as e increases. 
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FIG. 8. Conductance, as a function of n e , of a sin- 
gle-electron transistor with a = 0.25 and different values of e 
(T = 0.03£c). 



VII. DISCUSSION. 

In the following, we discuss some experimental evi- 
dence which can be explained within the model discussed 
here. 

It has been pointed out [Q that the correlations be- 
tween the conductances for neighboring charge states of a 



quantum dot |47[] are too weak to be explained using stan- 
dard methods for disordered, non-interacting systems. 
The experiments reported in [ p7|] are in the cotunneling 
regime. In the presence of non-equilibrium processes, we 
expect a behavior of the type g cx T 2 ~ 2e . Note that e is 
determined by phase shifts, which depend on microscopic 
details of the contacts. Thus, it can be expected to vary 
with the charge state of the dot, and to lead to large 
differences in the conductances of neighboring valleys. 

It has been shown that the temperature dependence of 
the conductance quantum dots, away from resonances, 
can be opposite to that expected in a system exhibit- 
ing Coulomb blockade fig ]. This effect could not be at- 
tributed to Kondo physics, as the data do not show an 
even-odd alternation. The reported behavior can be ex- 
plained within our model, assuming that the value of e is 
sufficiently large, and dependent on the charge state of 
the dot. 

In ref. 0] the inelastic contribution to the conduc- 
tance in a double dot sytem, where the electronic states 
in the two dots are separated by an energy e is measured. 
The result is approximately given by ii ne i(e) cx e A , where 
A is a negative constant of order unity. Taking into ac- 
count only one electronic state within each dot, the prob- 
lem can be reduced to that of a dissipative, biased two 
level system. The inelastic conductance reflects the na- 
ture of the low energy excitations coupled to the two 
level system. The observed power law decay implies that 
the spectral strength of the coupling, that is the function 
J(u>) in the standard literature |}43| , should be ohmic, 
J(ui) cx \u>\. This has led to the proposal that the exci- 
tations coupled to the charges in the double dot system 
are piezoelectric phonons |l!^,[50| . It is interesting to note 
that the excitation of electron-hole pairs leads also to an 
ohmic spectral function. Thus, at sufficiently low ener- 
gies, an orthogonality catastrophe due to electron-hole 
pairs shows a behavior indistinguishable from that aris- 
ing from piezoelectric phonons. The contributions from 
the two types of excitations can be distinguished at the 
natural cutoff scale for phonons, which is the energy of 
a phonon whose wavelength is of the order of the dimen- 
sions of the device. On the other hand, the simplest pre- 
diction for the expected behavior of the current induced 
by the emission of the electron-hole pairs is: 



1(e) cx K- 



1 



A 2 



1 



(31) 



where Ao is the tunneling clement between the two dots, 
et, is the bias, and K is the coupling constant (referred 
to as e in other sections of this paper). This expression 
gives the absorption rate of a dissipative two level system 
in the weak coupling regime, K <C 1 43 1. The natural 
cutoff for electron-hole pairs is bounded by the charging 
energy of the system. It would be interesting to disentan- 
gle the relative contributions of electron-hole pairs and 
piezoelectric phonons to the inelastic current. 

The photo- induced conductance in a double dot system 
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has also been measured |5T|] . The analysis of the contri- 
bution of inelastic processes due to electron hole pairs to 
the measured conductance proceeds in the same way as 
in the interpretation of the previous experiment. Let us 
suppose that a photon of energy Lu p h excites an electron 
within one dot. Assuming that the coupling to the envi- 
ronment is weak, the rate at which the electron tunnels to 
the second dot by losing an energy e is given by eq. ( 31 ) . 
The experiments in |5lj suggest that the number of states 
within each dot are discrete. Then, the induced conduc- 
tance should show a series of peaks, related to resonant 
photon absorption within one dot. The height of each 
peak is determined by the decay rate to lower excited 
states in the other dot, and it can be written as a sum 
of terms with the dependence given in eq.( |3l|) , where e is 
the energy difference between the initial and final states. 
The envelope of the spectrum should look like a power 
law, in qualitative agreement with the experiments. 

It is interesting to note that bunching of energy levels 
in quantum dots have been reported 153] . The separa- 
tion between peaks defines the charging energy, which, 
according to the experiments, vanishes for certain charge 
states. This behavior can be explained if the excitonic 
effects drive the quantum dot beyond the transition, and 
charging effects are totally suppressed. This mechanism 
can also play some role in the observed transitions in 
granular wires pq,p3| . 



VIII. CONCLUSIONS. 



We have analyzed the effects of non-equilibrium tran- 
sients after a tunneling process on the conductance of 
quantum dots. They are related to the change in the 
electrostatic potential of the dot upon the addition of a 
single electron. These effects can enhance or suppress 
the Coulomb blockade. The most striking effect arise 
from the formation of an exciton-like resonance at the 
Fermi level after the charging process and lead to the 
complete suppression of the Coulomb blockade and a di- 
verging conductance at low temperatures. It appears for 
sufficiently large values of the conductance, a , and the 
non-equilibrium phase shifts which define e in our model. 
The same potential which leads to this dynamic reso- 
nance plays a role in the deviations of the level spacings 
from the standard Coulomb blockade model [^2| . 

For simplicity, we have considered the simplest case, 
a single-electron transistor. The analysis reported here 
can be extended, in a straightforward fashion to other 
devices, like double quantum dots, where the effects de- 
scribed here should be easier to observe. 

As discussed, the non-equilibrium effects considered in 
this paper may have been observed already in the con- 
ductance of quantum dots. 
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